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ABSTRACT 


A model of the turbulent boundary layer separation 
of a two dimensional airfoil has been investigated. The model 
consists of aistributing sources on the contour in the separated 
region of sufficient strength to maintain the constant pressure 
characteristic of trailing edge boundary layer separation. A 
model for the attached boundary layer has also been included. 
An accurate estimation of maximum lift, pitching moment and 
pressure distribution have been obtained. A correct drag 
prediction has not been accomplished. 
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1. INTRODUCTION 


Analytical solutions to the problem of boundary layer 
separation have long been sought after but few have been success 
ful. The Navier-Stokes equations which govern flow separation 
have not been solved for high Reynolds number flows, hence 
existing analytical solutions must involve assumptions concern- 
ing the nature of separation. Empericism and/or experimental 
results are often used. In these cases, the question arises as to 
whether we are obtaining solutions which accurately reflect the 
physics of the problem, or are we modeling a phenomenon. This 
paper deals with the latter of these possibilities. 


One of the first and arguably the best analytical 
model of turbulent boundary layer separation was developed 
Jacob (Ref. 1). This model consists of using inviscid singe 
larity methods to model a viscous phenomenon. Sources are 
distributed in the dead air region to maintain constant pressure 
(in this region). A constant pressure, dead air region is a 
characteristic of turbulent boundary layer separation on an 
airfoil. The boundary layer and outer inviscid region must 
simultaneously be accounted for to complete the picture. Jacob 
has used this model to achieve very accurate prediction of 
maximum lift. 


A second model which has been used with similar success 
is the model proposed by Dvorak and Maskew (Ref. 2). Their model ' 
consists of placing vortices on the free shear layer. The 
strength of the vortices is such as to cbhtain constant pressure 
in the separated region. 


There are many other methods which have a stronger 
physical basis, but surprisingly, do not give results which can 


compare with the quality of results of Jacob and Dvorak's 
techniques. 


In this investigation, Jacob's model has heen chosen, 
because it has been in existance longer and has proven its 
validity as a model. 


It is the object of this investigation to develop the 
method to the point where an accurate estimation of the lift, 
drag, pitching moment and pressure distribution can be obtained 
for the turbulent boundary layer separation of a two dimensional 
airfoil. 


The importance of this investigation is that a correct 
prediction of the maximum lift coefficient is required. An 
inaccurate estimation of Ohaus can have drastic effects on the 
payload capacity of an aircraft (under estimation of Oo aie 
by 0.1 can mean up to a 20% reduction in payload). 


Two points must be emphasized before continuing. First, 
this method is only applicable as a design tool and will not 
give any insight into the mechanisms of separation. Secondly, 
separation is a three dimensional, and quite often unsteady 
phenomenon. This has been demonstrated by wind tunnel tests 
which have shown that under optimum conditions, it is difficult 
to obtain a truly two dimensional and steady measurement. Hence, 
the results of this research cannot directly be extended to a 
three dimensional situation. 
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THEORETICAL FOUNDATION OF THE MODEL 


2. 


1 The model 


The model is shown in figure 1. It consists of break- 
into three regions 


ing down the separated flow over an airfoil 


Region 1 : the outer inviscid region which is modelled by line 
vortices placed on the contour of the airfoil; 


Region 2 : the boundary layer characteristics which are deter- 
mined by an integral boundary method (similar to Le Foll's 

method). The boundary layer effect is included by distributing 
sources to compensate for the displacement thickness. 


Region 3 : the separated region is modelled by sources distri- 
buted on the contour of sufficient strength to maintain constant 
pressure on the airfoil in this region. 


2.2 The inviscid and separated regimes 
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The governing equation for the model can be derived 
from the boundary contours. Three boundary conditions are 
required for the analysis. Apptication of these conditions 
gives a Glauert type integral expression for the velocity 


on the contour 
(1) 
(2) 
(kinematic condition). 
(3) 
In addition to the boundary conditions, an additional criterium 
that the pressure is constant in the separated region must be 

enforced. Conditions 1&2 give rise to the integral expression 


for the vortex strength in function of blade geometry, uniform 
velocity, and the source distribution. 


Uniform stream far upstream. 


Normal velocity on the airfoil surface is zero 


An equivalent Kutta condition. 
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Martensen (Ref. 3) has shown that the flow over an airfoil may 
be considered as the superposition of four basic flows, a uni- 
form stream in the x-direction, a uniform stream in the 
y-direction, a circulatory flow, and an outflow. Expression 
(2.1) becomes four integral expressions for the vortex strength 
of the basic flows. 


Uniform stream in x-direction 


U 
Q 
< 


fa(t) [ry(zee) = -f (256) | dS(c¢) = & (2.2) 


Np |— 
a 


Uniform stream in y-direction 


§6(c) | (zie) & -#,(z.2) 2] aste) 
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Circulatory flow 


$u(é) | (2,0) 9% -¢ (240) 2 dS (x) 
y ds * dS 


Outflow : 
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Equation (2.1) has been split into one homogeneous and three 
non-homogeneous equations. The total solution is given by 


y(t) = Uly(c) - Vie(s) + Pu(c) + 8(c) (2.6) 
The following supplementary equations are required because the 


circulatory basic flow is the only basic flow which contributes 
to the circulation (note r in equation 2.6) 


$o(t)dS(t) = 0 (227) 
$8(c)dS(t) = 0 (2.8) 
gu(t)dS(c) = 1 (2.9) 
§6(c)dS(c) = 0 (2.10) 


Applying the equivalent Kutta condition determines the magni- 
tude of the total circulation : 


y(TEPS) = - y(TESS) (2511) 
Solving for the circulation : 


Uf a(TEPS)+a(TESS)} +V f 8(TEPS)+8(TESS)]} +{ 6(TEPS)+6(TESS)] 
C = + 
[ u(TEPS)+u(TESS)} 


The local velocity on the contour is given by 


W, = V¥(t)°#Q(c)? (2.99) 
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This expression contains a contribution due to the vortices and 
due to the sources (sources are distributed along the contour 
for the boundary layer model as well as the separation model). 


It appears that the complexity of the problem has been 
increased by separating the flow into four basic flows and that 
is true. The major advantage of this approach is that the velo- 
city on the contour for any angle of attack can be determined 
by substitution of the corresponding freestream velocity compo- 
nents in equation (2.6). One limitation is that if the source 
distributions is changed, the whole problem must be resolved. 


The solution of these integral equations is not 
straightforward. The equations become singular when calculating 
the influence of a particular vortex (or source) at its origin. 
This difficulty can be avoided as shown in appendix A. 


2.2.2 Numerical considerations 
Attention must now be turned to the numerical side 


of the problem. The profile is placed in the x-y plane as 
shown in figure 2. A coordinate transformation is used 


xX = (1-cos¢) (2.14) 
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This transformation allows a greater number of singularities 

to be placed near the leading and trailing edges where the 

slope of the contour changes rapidly. The contour is divided 
into N segments of constant included angle 4¢. On each segment, 
a line vortex of constant strength is placed. Line vortices 
have an advantage over point vortices in that a line vortex 


distribution more closely approximates the continuous vortex 


ae ae 


distribution of the integral equations. The integral expressions 
are replaced by the following summations 


Uniform stream in x-direction : 


1 N dy. 
— 2 a(t;5)K(z..05.85,) eee (2.15) 
Qn i=l J ds 
Uniform stream in y-direction 
1 N dx. 
— EF B(c,)K(z,.05,86,) = — (2.16) 
27 ji=1 ¥ dS 
Circulatory flow : 
1 N 
— u(o;)K(Z5505,40,) = 0 (23173 
2n i=l 
Outflow : 
LS ale, )K ) = LS ate.) ) (2.18) 
a BOLE IR(Z Ces hS.). Se SO Sa KZ eet shes 2. 
2n i=l ’ a A es Qn a ! ee a 
where the influence function K is given in appendix A. 
The supplementary conditions become 
N 
EB a(c4)aS(c5) = 0 (2.20) 
i=l 
N 
2 B(c;)dS(5;) = 0 (2.21) 
j=1 
N 
z u(o;)aS(z,) = ] (2.22) 


6(c;)aS(c;) = 0 (2.23) 
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Supposing that equations (2.15) through (2.19) must be satisfied 
in the points j=1,N where the vortices are located and together 
with equations (2.20) trrough (2.23) one obtains a set of 4(N+#1) 
equations with 4N unknowns. The unknowns are the vortex strength 
of the basic flows. However, a system of equations derived from 
an integral equation of the first kind is often ill-conditioned. 
This can be avoided if the kinetic condition is also applied at 
the intersection points of the line vortex. This leads to a set 
of 4(2N+1) equations with 4N unknowns, which can be treated by a 
least square technique. See Appendix A for further details. 


2.3 The boundary layer analysis 


2.3.1 The integral boundary layer approach 


The method used for the boundary layer region is by 
S. Huo (Ref. 4) and is similar to Le Foll's method. Huo's 
method is based on two assumptions 


(1) Local boundary layer velocity profile similarity. 


(2) The velocity profile is assumed to be a one parameter 
profile. 


The integral boundary layer equations are 


Momentum 
2 
2 * Me 
d(eu,) + 6 udu, = Ce a dS (2.24) 
Energy : 
u3 
ile = = cpuads (2.25) 


where the skin friction coefficient is 


2t 


Cr u2 
pie 


and the energy dissipation coefficient is 


c(y) 24 dy (2.27) 
ay 


The unknowns in equations (2.24) and (2.25) are 5°,0,0),C¢ 


and Cp: There are two equations at hand with five unknowns. 
In crder to solve these equations, the following information 
is necessary 


(1) A velocity profile law in the boundary layer. 
(2) A skin friction law. 
(3) A relation for the dissipation coefficient. 


The velocity profile law will enable the determination of the 
relation between the shape factors H,, and H35. Therefore, 
knowing two of the characteristic thicknesses of the boundary 
layer from equations (2.24) and (2.25) will allow the calculation 
of the other thicknesses. The skin friction and dissipation 
coefficients will be determined from the skin friction law and 
dissipation relation respectively. Hence, the complete boundary 
layer char-cteristics can be computed. 


Le Foil bases his method on the idea that the boundary 
layer can be described in terms of a form factor and a charac- 
teristic length instead of the velocity and length. The main 
advantage of this approach is that the inverse procedure of 
specifying boundary layer characteristics and their solving for 
the corresponding velocity distribution is easily achieved. The 
inverse problem is not of interest in this investigation, but 
the method is rapidly applicable to the direct problem. 


= 1053 


The parameters L (Truckernbrodt's parameter) and X 
are used by Le Foll to specify the boundary layer characteristics: 


dH3 2 


Eee Soe 
(Hi2-1) Hs2 


(2.28) 


As they are directly related to e-ch other, L and X can, there- 
fore, replace H,2 and Re, to characterize a boundary layer. 
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The assumed velocity profile is the Faulkner-Skan 
profile 


F(n,m) (2.29) 


Stewartson's transformation has been used with the Faulkner-Skan 
profile to extend the boundary layer calculation into the 
compressible regime. 


The compressibility aspect of boundary layer separa- 
tion is not included in the separation model, but it was felt 
that since a compressible boundary layer method was available 
and programmec), compressibility should be (and was) included 
in the boundary layer model. 


The skin friction and dissipation relation are func- 
tions of the assumed velocity profile. For the laminar boundary 


= IE 


To simplify the calculation, another change of 
variable is introduced 


Re ger 
foe a (2231) 
Re nef 
The momentum and energy integral equations become 
u C'(L) 
dP: S20 (Ee yds Se gS (2.32) 
ad rRe 
ref 
du 
GE. tS Bit ory 2 St (2.33) 
ds u. dS r dS 


The function C'(L,r) and B(L,r) are given in reference 5. 


Equations (2.32) and (2.33) are two first order dif- 
ferential equations which are solved by a fourth order Runge- 
Kutta procedure. 


2.3.3 Transition and separation 

Transition from laminar to turbulent flow can occur if 
the boundary layer becomes unstatle (dictated by Schlichting's 
stability curve). The unstable region is 


L + 0.011 - 0.005 X < 0 (2.34) 


Whether transition will actually occur depends on the turbulence 
level, pressure gradient, and Rey. A transition criterium has 
been developed by B. Roberts (Ref. 6) based on the experimental 
results of P.S. Granville. This criteria predicts transition if ARe, 
(the difference in Reynolds number between point in instability 
and point of transition) is larger than a specified amount 


(see Ref. 6). 
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The boundary layer may separate (laminar separation} 
before transition. Zero wall shear stress is used as the 
separation criteria. In the calculation of L, the integration 
constant is chosen so that L equals zero at separation. This 
corresponds to 


Hi> = 4.03 H39 = 1.515 


Since the inviscid model makes ro allowences for laminar sepa- 
ration regions, some approximations must be introduced. In all 
cases laminar separation bubbles followed by turbulent reattach- 
ment is suspected, the length of the bubble is neglected and a 
turbulent boundary layer is assumed from the point of laminar 
separation. 


2.3.4 The turbulent boundary layer 


The assumed velocity profile is a modified Falkner- 
Skan profile 


u 
ee ce g(n,2) (235)) 
Ue Us K 


The function g(n,Z) in equation (2.35) is determined from 
a combination of experimentai results and theory (see Ref. 5 for 
further details). The turbulent boundary layer is also appli- 
cable to compressible flow by the inclusion of Van Driest's 
generalized velocity. 


The skin friction coefficient is determined from 


ec 


T K 


_ (2.36) 


Ue a= S-an[ LPM? sncoante,. 
K 


a Pane 


The momentum and energy integral equations may be combined to 
obtain 


Co F(6 5086, (2.37) 


If or a given Mach number, similar solutions exist, a relation 
between the pressure gradient and velocity profile parameter G 
must also exist. Assume that 


2 
che Ee) ea ee (2.38) 
6.1 


where 


and 


|. 


u 
Ps, 
u { Hi2s 


with these assumptions, sufficient information is available to 
solve the turbulent boundary layer characteristics. See 
reference 5 for further details. 


Turbulent boundary layer separation is also indicated 
when L equals 0. At this point, the boundary layer calculation 
is stopped. 


2.4 The source distribution 


A source distribution is required for the separation 
region model. This distribution must be chosen in such a way 
as to maintain constant pressure (constant resultant velocity) 
in the separation region. Equation (2.1) which is in effect an 


expression for the resultant velocity on the airfoil indicates 


“a6 ae 


that the source distribution should be a function of : 
(1) profile geometry; 

(2) upstream condition; 

(3) angle of attack; 

(4) separation point location; 

(5) velocity distribution. 


This is a paradoxical situation. The source distribution depends 
on the velocity on the contour which in turn depends on the 
source distribution. Also, for a given airfoil, the fundamental 
relationship for the source distribution must be general enough 
to allow for changes in angle of attack. To change the source 
function for different angles of attack would defeat the purpose 
of using a functional relationship. 


Figure 3 shows a typical source distribution used in 
this study. The functional relationship is 


Q = 2.55"V 4, sina O< Ss" < 0.4 

Q= Veep. Sine 

Q = Veep, Sina- (Veg, Sina-0.2)(5S*-3)” / 4 

where 

5” is the distance along the contour in the separated region 
Veep. is the velocity in the separated region 


This relationship is dependent on angle of attack, velocity in 
the separated point and location of separation point. It does 
not include influence of the velocity distribution (over the 

remainder of the airfoil) on the separated region or influence 


of the profile geometry. A function which includes all of these 
aspects was too difficult to determine. 
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A completely different approach to the problem, which 
avoids functional source distributions altogether was attempted. 
This involves an iterative procedure to find a source distribu- 
tion maintained constant velocity in the separated region, but 
this procedure failed to converge. 


Another possible approach would be to a priori 
specify the velocity in the separated region, and then analyti- 
cally solve the inverse problem to find the source strength. 
There have been recent research efforts at VKI to model wing 
tunnel wall interference effects using source distributions 
determined by an inverse procedure. This research has shown 
that a solution of a nearly singular system of equations is 
necessary to find the appropriate source distribution. This 
sheds a Jittle light on the question why the iterative processes 
failed. 


In the unseparated region a source distribution for 
the boundary layer model is determined from 


Q = = (u.8*) (2.40) 


This relationship determines the source strength required to 
compensate for the displacement thickness. This approach (the 
blowing approach) is considered easier to use than the alterna- 
tive approach of using a displacement surface for a couple of 
reasons : the profile geometry does not have to be recalculated. 
The blowing approach is also more compatible with the separated 
flow model. 


2.5 Forces and moment on the airfoil 
2.5.1 The pressure coefficient 


The pressure coefficient is computed from the 


classical relationship 


Ww 
gc ie (2.41) 
p a2 
where 
P-P., 
ee pees q. = + pu? (2.42) 
Pog. 2 


These expressions are a result of Bernoulli's equation and are 
valid because this is an inviscid model. It is assumed that the 
pressure on the airfoil is the same as the pressure along the 
free streamline. Atso, no wake effects have been included. 


2.5.2 Lift coefficient 


The lift coefficient can be computed directly from 
the Kutta Joukowski law 


Pr c= $yds 
L = Teu_ (2.43) 
co 2r 


2.5.3 Drag coefficient 

An accurate determination of the drag for an airfoil 
with turbulent boundary layer separation is very difficult due 
to the foltowing 


(1) Integral boundary layer methods often do not predict the 
separation point exactly because these methods are not valid 
close to the separation points. 


(2) Singularity methods are not accurate near the airfoil 
trailing edge. 


(3) Surface roughness and freestream turbulence have a large 


influence. 
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(4) Wake effects have an effect especially when there is 
separation. 


(5) The equivalent Kutta condition does not model the physical 
Situation exactly. 


(6) Deviations from constant pressure in the separated region 
can have a significant effect. 


With these considerations in mind, the total drag is composed 
of a pressure contribution due tc incomplete pressure recovery, 
and a viscous contribution due to the viscous shearing stresses. 
The viscous drag is 


Qn 2 
D W 
Cy = y = cOSa | Ce(¢) as s d¢ 
Vv . puzc : de {ul 
2a ‘j 2 
+ sina | co(o) % [+] de (2.44) 
dod {u 
0 co 
The pressure drag is 
D on 
Cause —P__ = - cosa | C (¢) oy do 
p 1 pe dg 
2 
2a 
e BW ela ede (2.45) 
J Pp do 


The contribution of the separated region to the viscous drag is 
zero (Cr is zero in this region). 
2.5.4 Pitching moment coefficient 


The pitching moment is the integral of the pressure 


distribution over the contour and is given by 


EB: 


20 
ale dan a> | C(0.25-X) SX: dy 
4 

2 purc? 0 ae 

2 
en 

: | C_(0.25-x) 9X dg (2.46) 

0 : qe 


Expression (2.46) is the pitching moment coefficient about 
the quarter chord. 


2.6 The computer program 


Four computer programs have been written which perform 
the calculations described in the preceeding sections. The basic 
computational procedure is outlined in figure 4 and listings of 
the program along with examples of the input and output are in 
Appendix B. 


The calculations have been split into four programs 
to facilitate ease of use and also to allow modification of the 
results of each program if required. fuplication of computation 
is also avoided, for example, if for a given airfoil, the angle 
of attack is changed, the calculation may start with the second 
program (instead of the first program). 


3. RESULTS 


The separated turbulent boundary layer for the 
GA(W-1) and the NACA 2412 airfoil have been investigated. The 
results of these invesitgations are presented in figures 
5 through 18. 


For the GA(W-1) airfoil, figures 5 and 6 show the 
inviscid and viscid pressure coefficient for low angles of 
attack. Notice that the viscous calculated pressure coefficient 
agrees very well with the experimental results (Ref. 7) except 
near the leading edge wnere there is a short laminar separation 
bubble. 


At angles of attack greater than four degrees, compu- 
tational difficulties arose. The boundary layer program was 
originally written and calibrated for turbomachinery applications 
and in cases where the pressure peaks near the leading edge 
become too great, the boundary jayer calculation failed. A 
second problem was a short laminar separation bubble which 
existed throughout the angle of attack range. 


Figures 7-14 show the results of the calculation where 
the experimentel separation point has been used. The boundary 
layer model has not been included. The source distribution used 
for the separation region has heen optimized for an angle of 
attack of about 18 degrees, hence for other angles of attack, 
the pressure distribution is not as constant (in the separation 
region) as it is for 18 degrees. 


Figures 15 and 16 show the lift and pitching moment 
coefficient versus angle of attack respectively. Agreement with 
the experimental results is quite good. There is a deviation 
of the computed pitching moment coefficient at the high angles 
of attack. This is due to the calculated pressure distribution 
which differs slightly from the experimental pressure distribu- 
tion for the high angle of attack cases. Since the boundary 


layer program failed, a calculation of the drag was not achieved. 
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Figure 17 shows the results of a computed pressure 
distribution for NACA 2414 airfoil. A very flat pressure dis- 
tribution was obtained for this angle of attack. The lift 
curve for this airfoil is presented in figure 18. To obtain 
results for the NACA 2412 which included a prediction of the 
separation point the computational procedure was again modified. 


(1) The extremely high pressure peak near the leading edge was 
lowered (first five percent of airfoil). 


(2) An estimated separation point was used. 


(3) The boundary layer program was then used to validate the 
estimated separation point location. 


It was felt that lowering the pressure peak was justifiable 
because the boundary layer program predicted laminar separation 
in the first five per cent. A short laminar separation bubble 
would lower the pressure peak and since there is no laminar 


separation model, the peak was lowered by using engineering 
judgement. Figure 17 shows that a good prediction of Se 
was obtained. 


7 
4. CONCLUSIONS AND RECOMMENDATIONS 


The following can be concluded about the investigation 


(1) The model appears to be a valid model for turbulent boun- 
dary layer separation. 


(2) Could improve the method of determining the source distri- 
bution for the separated region. The inverse problem is a 
possibility. It will be very difficult to define an iterative 
procedure to find the source strength because of convergence 
problems. ‘ 


fi 
4 


(3) The boundary layer program should be extended so that it 
can handle airfoils with high pressure peaks and at high 
Reynolds number. 


(4) A laminar separation model should be included. 


(5) The lift, pitching moment, and pressure distribution have 
been accurately predicted for two airfoils. 


(6) The drag has not been predicted. 
(7) A wake model should be incorporated. 


(8) Extension of the computer model to compressible flow is 
worth considering. 
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APPENDIX A - DERIVATION AND SOLUTION 


OF THE INTEGRAL EQUATIONS 


Derivation of the integral equations 


This appendix contains a detailed derivation of 
integral equations. Attention is also given to the solution 


of these equations. 


The complex velocity at point x,y induced by a vortex 


and 


The complex velocity due to a distribution of vortices ona 


closed contour is 


Wiz) = - L ¢ 2) ascc) 
20 X+i¥ 


separating the complex velocity into its components 


viz) = - = § WE) as(c) (A.3) 
on X<+Y 


v(z) = 4 g Are} (A.4) 
en X“+Y 


A uniform stream must be superimposed 


iV 


@ 


oO se 


viz) =u, - Lg EE) ascc) (A.6) 
2 X°+Y 

itz) < Ag M2) as¢e) (A.7) 
a5 X aye 


Sources are placed on the contour for the boundary 
layer model. 


Wiz) = + Qe) (A.8) 
Qn X+iY 


For sources distributed on a contour, the complex 
velocity is 


tz) = tag SED as iey (A.9) 
Qn X44 


Superimposing the source distribution 


viz) = uj- Lg Wd as(eye Lg MQ) gs, (A.10) 
Qn X24¥2 Qn X24y2 

V(z) = Vt YS. as(zy- Leg TUE) ash (A.11) 
2n =X Ye an. XT #Y? 


The kinetic boundary condition that the velocity 
normal to the airfoil surface is zero is enforced. This is 
equivalent to stating that the surface is a streamline. 


aS ax dS ay dS 


6 Utz) ho viz) & (A.12) 
dS dS 


Introducing (A.10) and (A.11) one obtains 


Oe UY. Ss ee g Yx(2) as¢zy| oy 
ds [an X24 | ds 


QO. 
< 


+ ie § XQ(o) ds(c)| ey 
on NTH" { ds (A.13) 


~y gk. Poe ds(c)| ax 
X°+Y | 


tg wo] 


Since dS(z) is not a function of x,y the derivatives 
can be taken inside the integral. 


Te a a Ee ae f (2,0) 2% -# (2,0) a dS(z) 
Qn y ds * dS 


+ Fate) [Fy (ze) & -fy(2,0) B (A.14) 
2 2 ds d 


we OG ox 


Integral expression (A.14) can be split into one 
homogeneous and three non-homogeneous equations corresponding 
to the four basic flows described in chapter 2. 


1 § a(s) if (250) SS se esey 2h as(ey = & (A.15) 

2n y ds. > dS dS 

=f BCE) f (z.0) 2% -F (2,0) i dS(c) = 9% (A.16) 

2n y dS dS dS 

1 d d : 
oe $ u(c) [Fye8) Pe -f (2.5) a dS = 0 (K.17) 

1 dx _ dy 

a $ 6(&) Fy (z08) ae f (255) 7a dS(z) (A.18) 


= - + $ Qs) [Fx (eoe) ax ~F (253) a dS(z) 
2n ds dS 


The total vortex strength is the superposition of the four basic 
vortex strengths. 


y(t) = Uja(z) - Vi8(z) + Tu(s) + 6(c) (A.19) 


Equations (A.15) through (A.18) are obtained by substituting 
equation (A.19) into (A.14) and collecting like terms. 


The far field boundary condition (uniform stream far 
upstream) has been implicitly satisfied. This can be demonstra- 
ted by taking the limit of (A.10) and (A.11). As x approaches 
infinity, tae velocity components become equivalent to the 
imposed uniform stream. This result shows that the singularities 


do not have an influence far upstream. 
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The supplementary equations (as described in 
chapter 2) are 


ga(z)dS(z) = 0 (A.20) 
$8(s)dS(s) = 0 (A.21) 
gu(c)dS(t) = 1 (A.22) 
$8(c)dS(c) = 0 (A.23) 


The equivalent Kutta condition fixes the circulation 
y(TEPS) = - y(TESS) (A.24) 


Using equations (A.19 and A.24) gives 


-U[a(TEPS)+a(TESS)] +V_{ 8(TEPS)+8(TESS)} +{ 6(TEPS)+5(TESS)] 


[ u(TEPS)+u(TESS)] 
(A.25) 


Solution of the integral equations 


The solution of the integral equations is complicated 
by the fact that expressions (A.15) through (A.18) are singular 
integrals. Before tue discussion of the singular integrals, a 
coordinate transformation which allows greater numerical 
accuracy wil] be performed 


x = © (1-cose) (A.26) 
2 

Hence 

dy _ dy do dx _ dx do 


dS do dS ds dg dS 


/ 2 2 
dS(z) = 2 " | do(z) (A.27) 
y 


Similarity 


2 2 
dS(t) = / (= + | dg (zc) (A.28) 


A new variable will be introduced 


2 2 
y'(%) = y¥(&) Fl + | (A.29) 
3¢ ao 


This definition of the vortex strength has no physical signifi- 
cance, it is just used to simplify the derivation. This 
definition also applies to the basic flow vortex strengths. 


The singular integral 


Expressions (A.15) through (A.18) are singular integ- 
rals. This computational difficulty can be avoided by finding 
a function whose derivative is equal to the integral. Thus 
the integrals will become equivalent to the function evaluated 
between the limits of integration. The function I(z,c) must be 
determined such that 


dt (292). [' (2,2) 8% -#0¢2,0) 4 (A.30) 
dg (<) , dg * dg 


Make the following assumptions 
Ad A 


+ 


(1) y'(c) is constant over the distance 9- x to o+ 5 So that 
the vortex strength in (A.15) through (A.18) can be taken out 
from under the integral sign. 

aX ay 5 
(2) Also assume that Tae and IEE) are constant in that 


region. 


ii A cal 


- 29 - 


Define the following two functions 


Hence 


do(c) X24y2 

y _dn -y 9 
GF CSE) 3 doe) d>(c) 
do(c) Cas 


New multiply ag by C, and ay by C, where 


a een en Deen OY 


c, = dele) do(z dp(z) d¢(c) 
Ve aes d > 
ae 
dy(z) 
4s dx dy dn 


- 40(¢) do(z) do(z) de(s) 


Co = ; 
1 2 Zz 
do(c) d¢(z) 
After simplification, it can be shown that the function 
1(z,¢) in (A.30) has been obtained. 


aitzrt) 2 ¢, SF ye, 6 (A.31) 
do(c) do(s) d¢(z) 


1(z,t) = CyF(z,c) + C,G(z,<c) (A.32) 


0 2 


Hence the integral equations (A.15) through (A.18) reduce to 
the following summations 


1 N dy. 
— Z a'(c,)K(z,6;540,) = —9 (A.33) 
2a i=l dy 
1 N dx 
— £ B'(c;)K(z,6,548,) = a (A.34) 
2n i=l d¢ 
N 
= aCe Kizyeivhegyrs: 0 (A.35) 
20 i=l 
N 
1 ss Si (eK (eye op ne) (A.36) 
2n i=l 
ee Pee Oey We ey Chee EN aS (a) 
T d¢ ¢ 
where 
Ac; Ag; 
K(Z50; 5A.) = | Z50;+ ss -] 25637 lige (A.37) 
The supplementary conditions are 

N 

z a(c;)dS(z;) = 0 (A.38) 
j=l 

N 

£. B16,)4S(¢;) = 0 (A.39) 
j=1 

N 

E u(s,)aS(cs) = 1 (A.40) 
1= 

N 

z 8(¢,)aS(c,) = 0 (A.41) 
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Pressure Coefficient for GA(W-1) 
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Fig, 10, Pressure coefficient for GA(W-1) 
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Fig. ll. Pressure Coefficient for GA(W-1) 
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Fig, 12, Pressure Coefficient for GA(W-1) 
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Fig. 13, Pressure Coefficient for GA(W-1) 
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Fig, 15, Lift Coefficient for GA(W-1) 


Pitching Moment Coefficient for GA(W-1) 
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Fig. 18. Lift Coefficient for NACA 212 
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